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Abstract 



O 

The mass evolution of Schwarzschild black holes by the absorption of scalar fields is investigated 

5h 

in the scenario of the generalized Chaplygin gas (GCG). The GCG works as a unification picture 
of dark matter plus dark energy that naturally accelerates the expansion of the Universe. Through 

> 

elements of the quasi-stationary approach, we consider the mass evolution of Schwarzschild black 

r- 

holes accreted by non-minimally coupled cosmological scalar fields reproducing the dynamics of 
the GCG. As a scalar field non-minimally coupled to the metrics, such an exotic content has been 

O 

interconnected with accreting black holes. The black hole increasing masses by the absorption of 

the gas reflects some consistence of the accretion mechanism with the hypothesis of the primordial 

^ origin of supermassive black holes. Our results effectively show that the non-minimal coupling 

H 

with the GCG dark sector accelerates the increasing of black hole masses. Meanwhile some exotic 
features can also be depicted for specific ranges of the non-minimal coupling in which the GCG 
dynamics is substantially modified. 
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Accretion of matter by massive compact objects has been currently investigated in the 
background scenario of the general relativity [T]. Such accretion processes are candidates to 
mechanisms of formation of supermassive black holes (SMBH) in the center of most active 
galaxies [2j. In particular, it should follow some analogies with the process proposed by 
Salpeter and Zeldovich where galaxies and quasars could get some of their energy from 
processes of accretion [3]. 

Primordial black holes (PBH) are suggested as those created at the end of inflationary pe- 
riod. Among several physical processes which prescribe their production, one could suppose 
some kind of phase transition to the radiation era involving topological defects. Black holes 
with M < 10 15 g should not be observed today, as they may have evaporated by Hawking 
radiation [4] and contribute to the flux of cosmic rays. Those with M > 10 15 g could produce 
some observable effects at the present epoch [5] , in particular, the effect of gravitational mi- 
crolensing, or more generically, some significant contribution to the amount of dark matter 
[S]. When these PBH are formed with sufficiently large mass so that they do not evaporate 
completely, in case of matter and energy being absorbed in large quantities, these objects 
should be related to the supermassive black holes observed today [7j. 

These theories are all conceived in a scenario where observational evidences strongly 
suggest that the Universe is dominated by some dark energy component and evolves through 
a phase of accelerated expansion [8]. In such a context, the fine-tuning of the vacuum 
energy and the so-called coincidence problem are still open issues for which the standard 
AC DM model [S] does not provide a complete explanation. Although the quintessence 
models [TOl ITT] are in agreement with CMB data and provide a reasonable explanation 
to the coincidence problem, the corresponding equation of state, = p((p) / p{4>) , usually 
introduces some illogical constraints as, for instance, — 1 < < 0. Another class of dark 
energy model is resumed by the generalized Chaplygin gas (GCG) [T2T[1~4] . which is derived 
from the generalized model of D-branes. It synthesizes a unified dark sector through an 
exotic equation of state for dark matter and dark energy. The GCG is particularly relevant 
in this respect as it is recognized by observational constraints from CMB [15], supernova 
[T6HT8] , gravitational lensing surveys [19] , and gamma ray bursts [20J . Moreover, it has been 
shown that the GCG model can be accommodated within the usual formation mechanisms 
for large scale structures [121 El HZ] . 

Our simplifying hypothesis for the dark sector therefore sets a Universe completely filled 
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by the GCG underlying scalar field, which yields, through its gravitational impact, a viable 
unified scheme to mimic the presence of both dark energy and dark matter. 

The subjacent analysis considers the possibility of accretion of matter by massive black 
holes from minimal and non-minimal couplings to a real scalar field which drives the dy- 
namics of the GCG [T2l [T7j. The scalar field pervading all the Universe is assumed as an 
additional channel for black hole mass accretion. 

However, one should notice that the study of growing black hole masses in the presence of 
scalar field couplings to the metrics has been initiated before the discovery of the late time 
acceleration of the Universe and, thus, before the proposal of any dark energy accelerated 
expansion mechanism. Once inside the black hole horizon, this scalar field would be absorbed 
and produce the black hole mass increasing. The growing mass during the life of a PBH 
is something relevant and cannot be overlooked. We are mainly interested in modeling the 
black hole mass accretion by the GCG through the hypothesis of considering its non-minimal 
coupling to the metrics. 

In the above-mentioned context, the GCG model is characterized by an exotic equation 
of state p3J [15] given by, 

P=~— , (1) 
p a w 

which, for instance, can be deduced from a generalized Born-Infeld action [T3] . 

By solving the unperturbed equation for energy conservation in case of the equation of 
state introduced above, one obtains through a straightforward mathematical manipulations 



where po is the present energy density of the Universe, and A s = This model interpo- 



lates the dust dominate phase in the past, where p oc a -3 and the de-Sitter phase, p oc —p, at 
late times. This evolution is phenomenologically constrained by a and A s parameters, both 
positive, and with < a < 1. Obviously, a = corresponds to the ACDM model. The case 
a = 1 corresponds to the equation of state of the Chaplygin gas scenario [12] and is already 
ruled out by data [15] . Notice that for A s = 0, GCG always behaves as matter whereas 
for A s = 1, it behaves as a cosmo logical constant. Hence to use it as a unified candidate 
for dark matter and dark energy one has to exclude these two possibilities so that A s must 
lie in the range < A s < 1. In particular, for a flat universe, a = corresponds to the 
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(2) 
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ACDM model, and the best fit in agreement with the SNe-Ia datas happens for a = 0, 999 
and A = 0,79 (2]]. 

By assuming that the GCG fluid has a homogeneous time-dependent underlying scalar 
field [T2HHI ITT] , with an appropriate potential, one can define an action in the form of 

S = J d 4 x^{F(<j))R-d a <i)d a <j)-2V(<j))}. (3) 

where V ((f) is the self-interacting potential and F {(f)) features the coupling to the metrics. 
In the above action, units have been chosen so that 8nG — c — h — k — 1 and the metric 
signature adopted is ( — h ++)• 

By evaluating the above action, one obtains 

T ab = d a (j)d b (j) - 9 -f (<9 C 0«9 C + 2V) + V a V b F - g ab DF } (4) 

and the correlated Klein-Gordon equation described by 

□0 - V + \f'R = 0, (5) 

where the prime denotes derivation with respect to <fi. 

For an isotropic, homogeneous and spatially flat universe, one has 

ds 2 = -dt 2 + a 2 (t) (dx 2 + dy 2 + dz 2 ) , (6) 

and the temporal component of the Einstein equation gives the energy constraint through 
the following relation, 

3iJ(FiJ-F» = ^ + V(0), (7) 

where dot denotes differentiation with respect to the physical time. The modified Friedmann 
equation, for the space-space components of Eq. Q, can thus be written as 

— 2 ^F + ^F ,2 ^j H = 3(F + IF' 2 ) H 2 + 1 (1 + F") <p 2 -V-F' (V + H<p) , (8) 

where H = a/a is the Hubble rate. 
Minimal coupling 

By setting F(<p) = 1 one recovers the minimal coupling and the Klein-Gordon equation 
becomes 

+ 3H(f) + V ((f)) = 0. (9) 



Assuming that G ao = T ab , the energy-momentum tensor is then given by 

T ab = d a (j>d h 4> - 9 -f («9 C 0<9 C + 2V) , (10) 
from which the temporal and spatial components respectively gives the density and the 
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pressure as = \ + V {<j>), and p<p = \ — V (0). Using Eqs. (2) and (1 ), after some simple 
mathematical manipulations, one obtains the field equation for the GCG as 



cty(a) 1 / 1 (u] 
da fl y 4^o + i- 1 > 

Non-minimal coupling 

Models with F(<f>) = exhibit a plenty of singularities. Homogeneous and isotropic 
solutions passing from the F(<p) < to F((f>) > domains are indeed extremely unstable 
|22j . Thus, our study of the non- minimal coupling has been constrained by the solutions 
where F(<j>) > 0. 

In this case, the density and pressure are defined from the energy momentum tensor, 
Tab — (p + p)u a u b + pg a b, so that, p^ = —2H — 3H 2 , and p^ = 3H 2 , so that the temporal 
component of the Einstein equation is thus written as 

pF + 3H0F' = ^~ + V. (12) 



Assuming the simplifying notation with (f> = x and setting R = p — 3p, the Eq. (|5j) becomes 

(xx' + V')F' + 3HF' x- — (p-3p) = 0. (13) 



Inserting the redefined variables into Eqs. (|8j) and (12), and combining them with Eq. (13) 
results into 

(p + p)F = (1 + F")x 2 + F'x(x' -H). (14) 

Turning back to the previous notation, i. e. x = <p = aH^, one stays with 

3F F' Hn 

(P + P)^ = (1 + F"){<P\a)f + Yp <P\a)£ + FY(a) (15) 

and the potential can be found through Eq. ^ rewritten as 

F'V = P (F + 2F' 2 ) - (p + p) (f + - ^F> - V + \{1 + 2F")^(0'(a)) 2 , (16) 

where, from the beginning, we have assumed that the scalar field dependence on the scale 
factor a reproduces the background of the GCG. 
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FIG. 1: Dynamical dependence of the scalar field <fi on the scale parameter a, for the GCG 
model with fixed a = 0, 999 and A = 0, 79, in case of a non-minimal coupling parameterized 
by F(<f>) = 1 — 2£0 2 . Red lines indicate the £ < solutions, blue lines represent the £ > 
solutions, and the black line corresponds to the minimal coupling solution. 

Assuming a perturbative dependence of the coupling to the metrics on the scalar field 
through F((p) = 1 — 2£0 2 , where £ is the coupling constant between the scalar field <j) and 
the Ricci scalar, we have numerically solved the system of differential equations given by 



Eqs. (15) and (16) in terms of the scale factor, a. 

In Fig. [T]one can depict the scalar field dependence on a in case of — 1 < £ < 1. Notice 
that for the solutions where £ < (red lines), the scalar field qualitatively behaves as 
the minimally coupled solution for the GCG. Due to the non-linearity of the system of 
differential equations that we have solved, there is no continuous transition from negative to 
positive solutions for £. One shall notice the same qualitative behavior for values of £ > 1/3. 
However, for the interval where < £ < 1/3, which includes the conformal coupling, £ = 1/6, 
one has singularities when a tends to zero. 

As we have implicitly asked at the beginning of our exposition, does the picture of a 
Universe described by the GCG is consistent with the models for accretion of mass into 
black holes? In a black hole, the gravitational force is so intense that any form of matter or 
radiation is next to be drawn into it. In the region of the event horizon of a Schwarzschild 
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FIG. 2: The aspect of the potential running through the scale factor, in correspondence with 
Fig. [Tj Blue lines represent £ > solutions and red lines represent £ < solutions. 

black hole, the mass variation has two factors, the amount of matter and energy present 
at the border, which will be absorbed, and the evaporation by Hawking radiation. In spite 
of being studied for the case of ordinary matter, perfect fluid and scalar fields (23], until 
the discovery of accelerated expansion of the universe, the accretion of matter was not yet 
correlated to dark energy components. 

We have investigated the process of mass accretion and matter evaporation, through the 
above obtained results for the Chaplygin gas. In case of a scalar field in the vicinity of a black 
hole, assuming that on the horizon one has r = 2M, all the energy flux would be absorbed. 
According to the Bondi accretion [TJ, for which the field is stationary and spherically and 
symmetrically absorbed, the black hole mass increasing would be given by 



M 



dfir 2 T r = 47rr 2 T r , 



(17) 



r=2M 



The amount of matter/energy evaporated by Hawking radiation, in this case, is given by 





M 



M 2 '' 



(18) 



where /3 is a model dependent constant |24j . At first glance, and for simplicity, we have con- 
sidered non-rotating black holes in the Einstein frame. We have assumed the Schwarzschild 



metric as 



ds 7 



2M 



dr + 1 



2M 



-l 

dr 2 + r 2 (d9 2 + sin (9) 2 dip 2 ) 



(19) 



from which the spherically symmetric version of the Klein-Gordon equation obtained from 
Eq. ^ can be written as 



<9 2 1 

— - H 

dt 2 r 2 



1 - 



2M\ d 



dr 



2M 



Or 



1 - 



2M 



(20) 



At the vicinity of a black hole, the cosmological scalar field could generate a gravitational 
field much weaker than the gravitational field black hole. Although one can approximates the 
metrics by the Schwarzschild one in this region, a relation between the scalar field that was 
approximated from the infinity to the black hole horizon, namely the field absorbed by the 
black hole, is demanded. By assuming the quasi-stationary solution (25] for a slow-variation 
of the field in the nearby black hole, one should have 

'2Af ~ 



4> (r, t) « 0c 



v - r + 2M log 



(21) 



Substituting the above solution into Eq. (|20j), it is possible to extend the validity of the 
approximation to <p c ~ and V'((j) c ) ~ [26J. The field on the event horizon can thus be 
parameterized by the same cosmological field approximated from the infinity with some time 
delay, 

C (*) » 0oo + <M, (22) 

where 0oo and 0oo are constants obtained from the temporal solution of Eqs. (jsj) , (j7| and 
(l8j). The scale factor can be depicted from Eq. = p c h, so that 



c/iO 



^4, + 



a 3(l+a) 



v/3 



(23) 



Using this solution into Eq. Q one obtains 

2M ~ 



rpr 

1 n 



2M 1 + F 



(24) 



and the mass distribution in dependence on the scale factor becomes 

! F' 



M = 16vriir 



;i + f ) 



AM 



a. 
M 2 ' 



(25) 
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FIG. 3: Mass evolution of a black hole on the vicinity of the GCG. We have computed the 
effects for the coupling parameter in the interval —1 < £ < 1. The red lines show — 1 < £ < 
solutions and the blue ones show < £ < 1 solutions. The mass, M, is in units of Mo and 
the time, t, is in units of the Hubble time. 

For a minimal coupling, and considering the Hawking radiation, one thus get 

M = i6ttMVoo 2 - (26) 

Fig. [3] shows the mass distribution considering the Hawking radiation for a non-minimally 
coupled scalar field, using 1 < £ < —1 for M = 1, as a reference value. For 1/3 < £ < — 1 
the mass increases faster and indefinitely, but for £ > 1/3, one can notice some instability 
effect so that a decreasing mass behaviour is identified. 

The mass evolution and the scalar field dependence on the cosmological scale parameter, 
a, are semi-analytically obtained for a perturbative non-minimal coupling of the GCG un- 
derlying scalar field with the Ricci scalar of the metric. It reveals some qualitatively distinct 
behaviors for the black hole mass, M(t), and for the scalar field, <ft(a), when compared with 
the results for the minimal coupling case. For black hole masses smaller than a certain 
critical value assumed by the non-minimal coupling parameter, £ > 1/3, the accretion of the 
scalar field can lead to an exotic mass decreasing effect. 

By observing the correspondence between the solutions depicted from Figs. [T] and |3j 
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respectively for 0(a) and M(t), one can notice that the curves described by red lines for 
4>(a) follows the same qualitative behavior of the free GCG solution described by a black 
line. 

When the competition between dark matter and dark energy cosmological phases of the 
unifying picture of the GCG is dominated by dark matter, the black hole mass evolution 
reproduces the accelerated mass increasing (red lines) as usual for the non-minimal coupling 
involving ordinary matter fields. Some previous issues involving modified f(R) gravity 
models can also reproduce the same effects. Otherwise, the dark energy dominated phase 
of the GCG naturally accelerates the expansion of the Universe. For a specific range of the 
parameter £, the inclusion of the non-minimal coupling into the equations anticipates the 
dark energy phase of the scalar field as one can notice through the blue curves describing the 
dynamics of phi (a). For d(j)/da ~ at Fig. [TJ one has the anticipation of the cosmological 
constant phase, and for dcj)/da > one has a kind of phantom dark energy effect (with 
an exotically negative kinetic energy term that leads to the violation of the strong energy 
condition). The peculiar mass decreasing effect quoted above follows this kind of phantom 
dark energy dominated phase which takes the place of the suppressed dark matter dominated 
phase. Even if no phantom field is involved, the non-minimal coupling with the GCG unifying 
fluid provides such an effective picture for this. 

In more general lines, the mass evolution of black holes is governed by two simultaneously 
competing mechanisms. First of all, the Hawking radiation decreases the black hole mass 
due to the emission of thermal radiation. In parallel, the accretion of the surrounding 
available matter and energy usually increases the black hole mass. Possible evidences of 
primordial black holes at present shall depend on the balance of these processes and, as one 
could notice from our analysis, on the nature of the Universe's underlying scalar field, and 
on its couplings. The unexpected hypothesis of black hole mass decreasing because of the 
accretion of exotic (phantom) dark energy could change qualitatively the evolution of any 
black hole, implying, occasionally, into novel observational aspects for both astrophysical 
and PBH's. 

Reporting about the quasi-stationary approach, we have considered the mass evolution 
of Schwarzschild black holes in the presence of a non-minimally coupled cosmological scalar 
field that analytically resumes a unifying prescription for dark matter and dark energy of the 
Universe. Some imprints of such a unifying formulation onto the evolution of super-massive 
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black holes and, reciprocally, onto the cosmological behavior of the scalar field have been 
discussed. 

Such results summarize our investigation of some generic possibilities of non-minimal 
coupling of a dark sector scalar field to the Ricci scalar. As we have supposed, it could 
result into some changing on the dark sector dynamics, as well as, on novel ingredients to 
the mechanism of accretion of matter into black holes. We have found that for a non-minimal 
coupling drven by £ > some unexplained instabilities are present. It reinforces that only 
modified couplings parameterized by F((f>) > could be used. The analysis of the mass 
distribution of a black hole embedded by the GCG dynamical field allows for connecting 
PBH's to standard model black holes, in case of GCG absorption. Such result is in complete 
agreement with those found in the previous analysis [7J [26j EZj , where the mass of the black 
hole reproduces the increasing behavior predicted by the more fundamental models. In 
addition, it also reproduces some previous arguments in favor of the cosmology ruled by the 
GCG model. 
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